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We develop a formalism for treating coherent wave-packet dynamics of charge and spin carriers 
in degenerate and nearly degenerate bands. We consider the two-band case carefully in view of 
spintronics applications, where transitions between spin-split bands often occur even for relatively 
weak electromagnetic fields. We demonstrate that much of the semiclassical formalism developed 
for the single-band case can be generalized to multiple bands, and examine the nontrivial non- 
Abelian corrections arising from the additional degree of freedom. Along with the center of mass 
motion in crystal momentum and real space, one must also include a pseudo-spin to characterize 
the dynamics between the bands. We derive the wave packet energy up to the first order gradient 
correction and obtain the equations of motion for the real- and k-space center of the wave-packet, 
as well as for the pseudo-spin. These equations include the non-Abelian Berry curvature terms and 
a non-Abelian correction to the group velocity. As an example, we apply our formalism to describe 
coherent wave-packet evolution under the action of an electric field, demonstrating that it leads 
to electrical separation of up and down spins. A sizable separation will be observed, with a large 
degree of tunability, making this mechanism a practical method of generating a spin polarization. 
We then turn our attention to a magnetic field, where we recover Larmor precession, which cannot 
be obtained from a single-band point of view. In this case, the gradient energy correction can be 
regarded as due to a magnetic moment from the self-rotation of the wave-packet, and we calculate 
its value for the light holes in the spherical four-band Luttinger model. 

PACS numbers: 72.10.Bg, 72.20.My, 72.25.Dc, 73.23.-b 



I. INTRODUCTION 

It often happens, in transport phenomena, that one has 
to consider carrier dynamics in bands which are coupled 
together. This coupling arises either through strong in- 
terband scattering or as a result of the bands being degen- 
erate, or both. The nearly degenerate case is particularly 
relevant in transport theory as transitions often occur 
between bands even at relatively weak electromagnetic 
fields. Such situations include two-dimensional systems 
described by the Rashba Hamiltoniani with strong scat- 
tering, the doubly degenerate heavy and light hole bands 
in the Luttinger mode! 2 , which is frequently used to 
model the valence bands of bulk zincblende semiconduc- 
tors, and the conduction bands of wurtzite structures 3 . 
The case of nearly degenerate bands has not, to date, re- 
ceived the attention it deserves^, despite the important 
role played by such bands in semiconductor spintronics 
systems^, whether in dealing with spin currents 7 -, spin 
generation^ and relaxation^, or spin injection across a 
semiconductor interface^. 

Spintronics systems lend themselves to a semiclassi- 
cal treatment, as the external electromagnetic fields vary 
on scales that are considerably larger than atomic size. 
The semiclassical formalism has had much success in de- 
scribing carrier dynamics and transport phenomena in 
condensed matter physics. In the non-degenerate case, 
the carrier dynamics can be obtained semiclassically then 
combined with the Boltzmann equation to produce ac- 
curate descriptions of the transport properties of many 
materials. This approximation is used in the descriptions 
of cyclotron orbits, conduction in solids, the Hall effect 



and magnetoresistanceii. An essential application of the 
semiclassical model, which is specifically relevant to our 
discussion, is in treating external fields that are not rep- 
resented by bounded operators, so that a perturbative 
expansion will not converge^. The most common exam- 
ple is provided by uniform electric and magnetic fields, 
where the potential is linear in position. 

We therefore develop, in this paper, a semiclassical 
description of transport in degenerate and nearly de- 
generate bands. One of our main purposes is to ex- 
tend the semiclassical approach, as developed by Sun- 
daram and Niui 3 ., to the case of coupled Bloch bands, 
in order to take into account the spin degree of free- 
dom. We illustrate the underlying physics by treating 
two bands, without loss of generality. Two-band mod- 
els are frequently an adequate description of the conduc- 
tion bands of many semiconductors^. In experiments 
on spin transport in semiconductors the carriers have 
traditionally been electrons^, as the strong spin-orbit 
coupling in the valence band causes holes to lose spin 
information much faster™. However, in recent years re- 
search has also focused on spin currents in the valence 
bands of semiconductors 7 , with a degeneracy which is 
usually greater than two, and the formalism we outline 
is straightforwardly extended to multiple bands. 

To formulate a description of coherent transport in 
coupled bands we may no longer work with each band in- 
dividually but must instead treat the coupled-band man- 
ifold as a whole. The condition for our theory to be valid, 
which in the one-band case states that there must be no 
transitions out of that bandii, translates into the require- 
ment that there be no transitions out of the manifold un- 
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der consideration. We will consider a wave packet made 
up of two bands, which is a suitable description of coher- 
ent transport, when the density matrix has off-diagonal 
terms and the relative phase of the two wave functions 
plays a crucial role. This approach allows us to retain 
the notion of the real-space center of the wave-packet, 
r c , which remains well defined. Moreover, in extending 
the formalism to two bands we are able, in the presence 
of a magnetic field, to recover Larmor precession, which 
is not possible from a one band picture. The additional 
degree of freedom of the two-band system can be taken 
into account by defining a wave function with the Bloch 
periodicity in such a way as to incorporate both bands, 
which allows us to derive the dynamics from a single- 
band point of view. The coefficients of the bands can 
then be grouped into a vector which we shall call the 
pseudo-spin, the structure and dynamics of which makes 
clear the gauge structure of the problem. An interest- 
ing fact which will emerge from our analysis is that the 
effect of the external perturbations can be incorporated 
entirely into the Berry curvatures^, which in turn are 
generated by a set of connections in real and reciprocal 
space as well as in time. The Berry curvatures acquire 
additional terms needed to ensure gauge covariance, and 
in the framework we present they take the form of field 
strength tensors associated with the connections. 

The organization of this paper is as follows. In Section 

II we develop the semiclassical formalism for coherent 
transport in the presence of electromagnetic fields, de- 
riving the Lagrangian, based on a time-dependent varia- 
tional principle, and the equations of motion. In Section 

III we use our formalism to show how coherent wave- 
packet evolution under the action of an electric field leads 
to the separation of up and down spins. This idea is sim- 
ilar in principle to the spin transistor proposed by Datta 
and Dasi&. We demonstrate that a large degree of tun- 
ability can be achieved by varying the gate field and num- 
ber density. Finally, in section IV we examine the case of 
a magnetic field. We show that the gradient correction 
to the energy can be interpreted as an intrinsic magnetic 
moment of the wave-packe1ii2iiLi& ! and we calculate this 
angular momentum correction for the light holes in the 
spherical four-band model of the Luttinger Hamiltonian. 



II. FORMALISM 

The semiclassical model describes the dynamics of 
wave packets. The wave packet we consider is well lo- 
calized in reciprocal space, and it is assumed it sees only 
a small part of the lattice at any one time. It is chosen in 
such a way that its spread in wave vector is much smaller 
than the size of the Brillouin zone, so that its motion at 
any moment is dependent only on the local properties of 
the band structure. In order for this to happen, the un- 
certainty principle dictates that the spread in real space 
must be greater than the size of the lattice constant. 

We consider systems whose Hamiltonians are func- 



tions of slowly varying parameters, such as the poten- 
tials of weak external electromagnetic fields, which vary 
on larger length scales than that of the wave packet, and 
are treated classically. The periodic potential of the ions 
on the other hand, changing over dimensions small com- 
pared to the wave packet spread, must be treated quan- 
tum mechanicallyii. 

Given these conditions, we define the local Hamilto- 
nian H c (r c ,t) as the Hamiltonian with the slowly vary- 
ing potentials evaluated at the center of the wave packet, 
which we denote by r c , and time t. The Hamiltonian 
may be expanded^ about r c and if the external fields 
vary on spatial scales much larger than that of the wave 
packet we may truncate the expansion at the gradient 
term, which we define by AH: 

1 riff 
AH=-[(r-r c ).^+c.c.}. (1) 

The gradient term gives rise to a correction to the energy, 
which will play an important role in our discussion below. 

The energy spectrum of the local Hamiltonian H c con- 
sists, as usual, of a series of bands, of which several are 
close together in energy and are separated from the others 
by larger gaps. It is the subset spanned by these bands 
that constitutes the focus of our attention. We regard the 
fields in this problem as small enough that Zener tunnel- 
ing to the remote bands is negligible, but they may still 
be strong enough to induce transitions within the sub- 
set. For an energy spectrum with such a structure we 
may further decompose the local Hamiltonian into a de- 
generate part, Hd, which, when restricted to the subset 
of bands closely spaced in energy is proportional to the 
identity matrix, and a non-degenerate part, H n , which 
is assumed small and treated perturbatively. The local 
Hamiltonian of is then: 

H c = H d + H n . (2) 

The gradient correction to H c can also be expressed in 
terms of the degenerate and non-degenerate contribu- 
tions: 

AH c = AH d + AH n . (3) 

Since H n is treated as a perturbation the gradient cor- 
rection to it, AH n will be second order in smallness. We 
will therefore neglect this correction henceforth. 

When the external fields are smoothly varying the 
states move within the subset of bands which are close in 
energy and which henceforth, for simplicity and with- 
out loss of generality, we take to be two-dimensional. 
The subset is spanned by two basis functions, which are 
eigenstates of Hd , the degenerate part of the local Hamil- 
tonian, evaluated at r c , which has the periodicity of the 
unperturbed crystal: 

H d |¥ i (r c ,q,t)) = e|* i (r c ,q J t)). (4) 
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For a given r c , therefore, these eigenstates have the Bloch 
form, with the functions |uj) representing the lattice pe- 
riodic parts of the wave functions: 



|*i(r CJ <Lt)>=e*»-'|«i(r CJ cbt)) 
|* 2 (r c ,q,t)) = e lq > 2 (r c ,q,i)) 



(5) 
(0) 



The wave functions |ui(r c ,q, t)) are spinors with the full 
periodicity of the lattice. Despite the fact that the two 
bands are spin split, it cannot be assumed that their lo- 
cal spin quantization axes are antiparallel, as the interac- 
tions with neighboring bands may affect the direction of 
quantization. Therefore, in principle, a finite overlap ex- 
ists between the spinors corresponding to the two bands 
and it is not revealing to make a further decomposition 
of the eigenfunctions into an orbital and a spin part. Ad- 
ditionally, the Hamiltonian contains terms describing the 
spin-orbit interaction, which may depend on wave vector 
and position. 

Employing the crystal momentum representation, the 
wave packet is expanded in the basis of Bloch eigenstates: 



= J A{a(q,<)[^i(q t)l*i>+»h(q.*)l*2)]}. 



(7) 



As the wave packet depends only on the local proper- 
ties of the band structure, the basis functions l^i), l^) 
are functions of the position of the wave packet cen- 
ter, r c , wave vector and time, although implicit in the 
ket notation is dependence on position. The function 
a(q, t) = |a(q, t)|e~ ir ( q,t )' 2 , which incorporates the over- 
all phase term, is a narrow distribution function describ- 
ing the extent of the wave packet in reciprocal space and 
is sharply peaked at the center of the wave packet, de- 
noted by q c , as discussed by Sundaram and Niu 13 . The 
functions 771 and 772 describe the composition of the wave 
packet in terms of the two bands. The wave packet sat- 
isfies the normalization conditions: 



A|a| 2 = l,|*?i| 



= 1. 



(8) 



The wave packet can be rewritten by grouping together 
the coefficients in an overall wave function |u), which 
retains the Bloch periodicity: 



<Fq\a\e 



ir/ V q - f |«). 



(9) 



Note that |u) is not an eigenstate of the local Hamiltonian 
H c = Hd + H n , but an expansion in eigenstates of Hd, a 
crucial difference from the one-band situation. In addi- 
tion, the wave vector and time dependence of |u) come 
both from the time dependence of the Bloch states and 
that of the coefficients. 

We require the real space center of the wave packet to 
be given by: 



{w r \w) 



dT c 
<9q c 



(10) 



The subscript c signifies that the quantity is evaluated 
at the center of the wave packet in reciprocal space, that 
is q = q c . The vector R, representing a connection in 
reciprocal space, is defined as follows: 



d 

R = ^dq'"^ 



(11) 



The energy of the wave packet is given by the expec- 
tation value 



(w\H\w) = (w\H d \w) + (w\H n \w) + (w\AH d \w) 

= 6 + A n + A d = £. 



(12) 



Both A n and Ad are expressible entirely in terms of the 
Bloch wave function |u). A n is given by 



a„ = (u\H n \u) = vlKm, Km = 



while Ad is 



1 , , , dHd du 



dr c <9q 



c. c .)--.R. 



(13) 



(14) 



In the above H n = e _iq, *H n e iq -* while H d = e- iqf H d e iqf . 
The energy correction Ad is identical to the expression 
obtained by Sundaram and Niui^. It takes on an addi- 
tional significance when a magnetic field is present, as 
will be seen in the last section. 

The Lagrangian C is obtained semiclassically by means 
of a variational principle: 



C = (w\{ih--H)\w) 



(15) 



Its use is justified by the fact that the Euler-Lagrange 
equation of motion for |w) derived from it is the time- 
dependent Schrodinger equation. Following the method 
used by Sundaram and Niu, the following expression is 
found for the Lagrangian: 



C = (u\ih—) + hr, 



q c - e = 

hx c ■ q c - 
77* (A" +Af J ) Vj . 



(16) 



In the above 4^ represents the total time derivative, in- 
cluding both the explicit time dependence and the im- 
plicit, which is due to dependence on q c and r c . The La- 
grangian depends only on the values of rji and along 
the trajectory q = q c (t). Since q c is a function of time 
only, we may regard rji in the Lagrangian as an inde- 
pendent variable, f)i(t). The equations of motion derived 
from the Lagrangian are: 



ftq c = -q^- + (^r r r c + fJ rq q c ) - flu 

de 

hv c = — - (f2 qr r c + fl qq q c ) + fi tq 



(17) 
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Here Hij = (uj|H c |uj). The curvature tensor fi™/ is de- 
fined by: 



du . du 



du , du 



x dr a drp 
and the vector tt t q by: 

, du . du . 



tq 



du . 9u 



dt dq a dq a dt 



)) 



(18) 



(19) 



The others can be deduced analogously. These quantities 
have exactly the same form as the curvatures defined in 
the paper by Sundaram and NiuA£. 

We specialize in the case of an external electro- 
magnetic field. The effect of such an external field is 
discussed thoroughly by Sundaram and Niuii The wave 
vector q must be replaced by k = q + f A(r, t), which 
is the gauge invariant crystal momentum (for electrons 
with charge -e), and therefore the Hamiltonian will 
have the form H(k) +eV(r,t). Provided the magnetic 
or exchange field is constant and uniform, so that the 
Zeeman term has no time or space dependence, the 
basis states {|u;)} will depend only on k. The reason 
for this is that all the spatial and time dependence 
of the wave functions will only come from the spatial 
and time dependence of the vector potential A(r,t). 
We will therefore restrict our attention to constant 
uniform magnetic fields, while the electric fields may 
be space- and time-dependent. As the electromag- 
netic fields vary on a spatial scale which is large 
compared to that of the wave packet, the local Hamil- 
tonian will have the form H[q + ^A(r c , t)] + eV(r c , t). 
The band eigenstates {l^nk)} take the form 



IV'nk 



3 iqr 



U„k) 



Unk; 



The time de- 



pendence of |u) comes both from the Bloch wave 
functions {|ui)}, which depend only on k, and from the 
coefficients, which depend only on time. Therefore, the 
Lagrangian in the presence of electromagnetic fields can 
be written as: 



C = ^^dt^ + { hkc ~ eA ( rc ' *)] ' / 2Q j 
-e-A n -A d -eV(r c ,t). 

The equations of motion now take the following form: 



dk r 



hk c = -e(E + f c x B) 
(u\H\u) - Hk c x $7 + fi tk 



= (Hij - ftk c (uj|i^-»%, 

where fi = i(§j| X |f^)- Note that the position vector 
equation of motion is very similar to the one band casein, 
excepting the presence of the vector f2 t k, which is non- 
zero due to the time dependence of |u) through the co- 
efficients. The equation of motion for |u), if a magnetic 
field is present, leads to the formula for Larmor preces- 
sion. The equations may be solved to any desired order 



in the external fields and are not limited to the linear 
response regime (the fields are weak enough that they do 
not induce transitions to remote bands). 



III. THE PROBABILITY AMPLITUDES 

The treatment we have presented so far is an exact 
analogy with the single-band dynamics. The equations 
of motion l|21|l are complete. Nevertheless, the equations 
of motion can be made more explicit in terms of the co- 
efficients ?7i, and the non-Abelian quantities emerging in 
the process illustrate the gauge structure of the Hilbcrt 
space. 

The coefficients 771 , 772 give the composition of the wave 
packet in terms of the two bands, and it is natural to 

think of them as a vector, ( ^ 1 ) , which will be called 77. 

The connection R can be expanded in terms of rj: 



t 9i] 



, where JZ? 4 



(uilip-)- (21) 



We will also introduce the time connection 7y = (uj |i^-) ■ 
The Lagrangian in this picture takes the form: 



C = ihrj 



151 



r^TL-q. 



(22) 



where the covariant derivative with respect to time, de- 
fined as j^r = ^ — i(T + q c • R), has been introduced. 
Specializing in electromagnetic fields, we end up with the 
following Lagrangian: 

C = rf(ih—)r) + [hk c - eA(r c , t)] • r c - 
-tfHr]-eV(r c ,t). 

The equations of motion derived from this electromag- 
netic Lagrangian are as follows: 



D 



-e(E + r c x B) 



(24) 



r7 t [— — , Ti]ri - hk c x r/^Ti] 
Dk 

■fi Dr > v 
in—— = nri 

Dt 1 

The covariant derivative with respect to the wave 
vector has been introduced, which has the form 
= g| iJZ a . The non-Abelian Berry curvature ma- 
trix, J-?, is expressed in terms of the field strength tensor 
corresponding to the covariant wave vector derivatives: 



where 



T 



dm 



T 1 - 



a/3 



1 



a/37 j: a @ 



D D 



l Dk a ' Dkp 3 



(25) 



(26) 



dn° 



dk a dkp 
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This form, which includes the non-Abelian correction 
from the commutator of the connection matrices, makes 
evident its gauge covariance with respect to unitary 
transformations of 77. The curvature tensor is antisym- 
metric under interchange of a and /3, while the indices i 
and j satisfy = {Jy^)* . 

It is seen from the equations of motion that working in 
the coupled-band manifold entails the presence of non- 
Abelian quantities such as the modified Berry curvature 
and gauge covariant group velocity which are 

corrections to the one band equations of motion needed to 
ensure gauge covariance. The matrix TL is not necessarily 
diagonal, as it may include energy gradient corrections. 

We note that equivalent results can be derived using 
an argument based on the Ehrenfest theorem, as in the 
extensive work of Shindou and Imura 22 . 



IV. CONSTANT ELECTRIC FIELD 

We will examine first the case of a constant uniform 
electric field acting on two degenerate bands. We choose 
a gauge such that the scalar electric potential need not 
be included in the Hamiltonian, and the electric field 
is represented purely by the vector potential A. With 
experiment in mind, we take E=(0,0,E), modeling a gate 
field, and study its effect on transport in the xy-planc. 



A. Electrical spin separation 

We choose as an example the spherical four-band 
model: 

H Lutt = |1[( 71 + L 2 )k 2 - 2 72 (k • J) 2 ], (27) 
2m 2 

where J is the total angular momentum operator, m 
is the bare electron mass and 71 and 72 are material- 
specific parameters. The wave functions are eigen- 
states of the helicity operator k • J and have the form 
|u m ) = c _1 * Jz e _l6,Jy |m) where |m) are eigenstates of the 
angular momentum operator J z while 9 and <f> are the 
polar and azimuthal angles of the wave vector, respec- 
tively. We shall treat the two-fold degenerate heavy and 
light hole manifolds separately and we shall denote the 
probability amplitudes in the heavy hole subspace by 77 11 
and those in the light hole subspace by rj L . 

In these subspaces, the equations of motion for the 
probability amplitudes take the form 



IT 



if&l = (H H - eETZ zH ^ H 

.,d V L 
in—— 

dt 



(28) 



{H h - eEn zh )ri h , 



where the superscripts H and L represent restrictions to 
the heavy and light hole subspaces, respectively. The re- 
ciprocal space connection matrix is given by the following 



expression: 



11 



09 



J y 



-{J z cos( 



J x smt 



(29) 



Ok Ok' 

In the heavy hole sector TZ z =0 and the bands decouple, 
therefore no spin separation can be achieved electrically 
in the heavy hole manifold. Henceforth we shall concen- 
trate on the light hole manifold, where the connection 
matrix 1Z Z = — ^r& y has off-diagonal elements only, with 
k^ = (k x ,k y ) and er y a Pauli spin matrix. We shall sup- 
press the index L in what follows. 

The equations of motion for the position and wave vec- 
tor are: 

hk = eE 

dei „ + _ 
hr c = — - eE x rj'Ti], 
ok 

in which ko is the initial value of k, £1 — 4-^- is the light 
hole energy, mi is the light hole effective mass, and the 
curvature T = f p-c z - The wave vector equation of mo- 
tion is readily integrated to give k = ko + . Since the 
Berry curvature vector is parallel to k, there are two lim- 
iting cases to consider: the case k //E is trivial because 
the curvature correction vanishes and the bands decou- 
ple, so we will focus on the more interesting case ko _L E. 

The equations of motion can be solved exactly, rj is 
given by: 



/ (o) , (o) . > 
_ / 1][ cos a + rj 2 sin a 

V - (o) (0) . 

\ 772 cos a — r/l sin a . 



(30) 



with the angle a(r) = arctaji( T +™ s o 9 ° ) - (§ - 6 ), where 
we have introduced the dimensionless time r = §1^ and 

9q is the polar angle of ko , and where r]> are the values 
of rj at r=0. 

In this system, the contraction ry a Y r\ (with i=l,2,3) is 
the expectation value of the pseudo-spin. Its components 
evolve in time as: 

(a 1 ) = (<r 1 ) T = cos2a — (a 3 ) T= o sin 2a 

(^ 2 ) = (a 2 )r=o 

(a 3 ) — (cr 3 ) r= o cos 2a + (<t 1 ) t= o sin 2a. 

The electric field therefore only rotates the 1 and 3 com- 
ponents of the pseudo-spin into combinations of each 
other, while the 2 component remains unaffected. To 
understand the significance of these results we will ex- 
amine a concrete example, taking initially a positive he- 



licity eigenstate so that = 1, 772"'' = 0, and fixing the 
initial wave vector along the x-axis such that ko = kox, 
which means that 9o — ^ . The full time evolution of the 
pseudo-spin components is: 



(0) 



1-T 2 



2r 



{a ) = {a ) T=0: 



- {a- J ) T =o-_ 



(a 2 ) = (<7 2 ) T=0 (31) 



2t 



\P )t=0: 
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As t — > oo, a reaches the limiting value of ^ and the 
components of the pseudo-spin become: 



{<J ) T =0 
(<T 2 ) r=0 
(<T 3 }r=0. 



(32) 



Thus the 1 and 3 components of the pseudo-spin are re- 
versed while the 2 component is conserved. 

The time evolution of the wave vector is described en- 
tirely in terms of the time evolution of the angle 9, which 
is most conveniently expressed as 



cos 9 = — 
k 



sin 9 = — — = 

k 



< El 

h 



^o 2 + ( £ f) 2 ^ 
k Q 



72' 



(33) 



Therefore, initially we have cos 9=0 and sin 0=1 while as 
t — ► oo , cos 9 — >1 and sin9 — *0. 

The expectation value of a spin component operator 
s Q in the wave packet |w) is given by (w|s Q |w) = rj's a T], 
where Sy = (ui|s Q |uj). The time evolution of the spin 
of one electron can thus be found by knowing the time 
evolution of its pseudo-spin. Since our goal is to sepa- 
rate spins of opposite orientations, it is sufficient to know 
the value of the pseudo-spin. The bands being spin-split 
holes with pseudo-spin up also have spin up and holes 
with pseudo-spin down have spin down. However, it is 
instructive to follow the motion of the spin as time pro- 
gresses, as well as the time evolution of the helicity. The 
expectation values of s x , s y and s z are: 

(s x ) = — (— sin 9 cos 0(<t 3 ) + cos 9 cos ^((T 1 ) — singer 2 ) 

(s v ) = — (- sin6'sin0((T 3 ) + cos 9 sin cf> (a 1 ) + cosher 2 ) 
3 2 

(^) = ^(icos^ 3 )-sin^ 1 )). 

(34) 

We assume the carriers have been polarized (by opti- 
cal means, for example as done in the experiments of 
Malajovich et aZ i 19 i 20 , although those utilized electrons) 

so that r/ ' is either or Kj^ . Therefore the initial 

expectation values of {o -1 } are: 



(35) 



and the initial spin expectation values are given by 



(a 1 ; 


>T = 




= 




> T = 




= 


- 3 >T 




(- 3 >4 


= 1, 



(J*)o = ±-,(^) = (s z }o = 0. 
o 



(36) 



It can then easily be seen that the y-component of the 
spin is zero at all times. Substituting for 9 and {a 1 } in 
B34JI we obtain the time evolution of the other two spin 
components: 



<^>0 = ±7T 



1-5t 2 



3 2(l + r 2 ) 3 /2 



3 2(1 



)3/2 



(37) 



As r — > oo, the expectation values of the spin components 
are: 



(38) 



The spin in this case is not conserved. However, a closer 
look at (|34|l reveals that (s x ), (s y ) and (s z ) cannot be ob- 
tained from (a 1 ), (<r 2 ) and (a 3 ) by a rotation, as the ma- 
trix describing the transformation is not unitary. There- 
fore one should not think of the projection of the spin 
onto the light-hole subspace as a vector. 
Finally, the helicity is given by: 



k-(s) 



-(a 



hi 
± — 

6 1 



(39) 



The helicity is proportional to the expectation value of 
the third component of the pseudo-spin. It is therefore 
not conserved for light holes in an electric field. This 
conclusion was also reached by Jiang et a/iSi. 

The r c equation of motion can be integrated to give 
the trajectories of the carriers: 



ft 2 fcg , „ 

— (rx - 



eEmi 



r(3 + r 2 )(| ?? (")| 2 ^| 7? (°)| 2) 
2fc (l + r 2 ) 3 / 2 



(40) 

We have omitted a term proportional to il^r]^ since in 
our setup either one of them will be zero . The second 
term in l|40|l will have opposite signs for the carriers with 
rj initially up and those with r\ initially down. Therefore, 
these carriers will be separated in the y-direction. From 
the above and Fig. 4.1 it can be seen that the maximum 
separation in the y-direction occurs at r = 1 while as 
t — > oo this separation tends to l/k . 



B. Experimental observation 

We discuss an experimental setup in which the effect 
we have described can be measured. We propose using a 
three-dimensional semiconductor slab containing a non- 
degenerate hole gas. The sample must be clean in order 
for the hole spin relaxation time to be long, specifically 
of the order of picoseconds. Carriers are excited optically 
from the conduction band into the valence bands by using 
a laser beam. Provided the laser beam is sharp, only a 
narrow range of k-space will be excited around k=0. The 
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optically excited holes will have wave vectors lying in a 
narrow spot about the origin. We assume they have been 
excited into a state of definite spin. Both light and heavy 
holes are excited but, as shown in the previous section, 
the heavy holes do not separate according to spin under 
the action of the electric field. A source and a drain will 
be positioned along the x-direction on the two faces of 
the sample while a gate terminal will be present on top. 
After the optical excitation, the magnitude of the holes' 
wave vector can be increased by applying a source-drain 
field E x in the form of a picosecond pulse, which will ac- 
celerate the carriers along the x-axis, its magnitude tuned 
to ensure ko has the desired value. This source-drain field 
provides an additional advantage. In the process of opti- 
cal excitation electrons as well as holes will be excited in 
the sample and the field which drives the holes one way 
will drive the electrons the other way, ensuring the effect 
observed is indeed due to holes. By adjusting the mag- 
nitude of the source-drain electric field pulse the initial 
wave vector ko of the holes incident upon the interface is 
tunable over several orders of magnitude. We will choose 
a source-drain electric field in such a way that the wave 
vector ko will have an x-component which overwhelms 
the y- and z-components. We will also choose the mag- 
nitude of ko to be approximately 1/b, where b is the real 
space thickness of the laser beam. The reason for this is 
that in the limit of large r the spins are separated by a 
distance of approximately 1/ko, therefore the separation 
of the spins will be approximately the same as the width 
of the laser beam. Once excited, the carriers will be sub- 
jected to the action of the gate electric field E along z, 
which will lead to separation of spins as described above. 
The spin accumulation at the other end of the sample 
can be measured by Faraday or Kerr rotation. It will be 
position-dependent along y, that is, as one moves along 
y the spin-z polarization will change sign. 

We take the dimensions of the slab to to be 50 nm x 
5 /j,m x 5 /jm and the width of the laser beam is taken 
as 1 /im. The optically excited holes will be accelerated 
until their wave vector reaches the value fco = 10 6 m . 
For a source-drain field E x of 500 Vm _1 and a light hole 
mass of O.lmo, where mo is the bare electron mass, the 
distance traveled by the light holes along the x-axis will 



be 



2micE x 



7.2 nm. This will happen after a time of 



1.25 ps, which can be achieved in samples in which the 
holes have longer spin lifetimes. Therefore the source- 
drain field must be a 1.25 ps pulse of amplitude 500 V 
m . 

We will take the gate electric field E=25000 V m _x . If 
one waits for the value of r to reach 50, then the mag- 
nitude of the spin polarization along the z-direction will 
be approximately h/6 while along the x-direction it will 
be negligible. The separation between the carriers with 
spin-z up and spin-z down will be approximately 1 /xm, 
which is the same as the real-space width of the laser 
beam and thus observable. The waiting time will be ap- 
proximately 1.3 ps. Finally, the distances traveled in the 
x- and z-directions under the action of the gate electric 



field are 35 nm and 850 nm respectively. 

This phenomenon is similar to effects such as the spin 
Hall effect since carriers with different helicities are sep- 
arated in the xy-plane by the electric field normal to the 
plane. 




1 2 3 4 5 

Dimensionless time 

FIG. 1: Separation 1 in the y-direction between light holes of 
opposite helicities as a function of r, the dimensionless time. 



V. CONSTANT MAGNETIC FIELD 

When a constant uniform magnetic field is present, the 
gradient correction to the degenerate part of the Hamil- 
tonian takes the form: 



A rf = — M • B. 



(41) 



M, which is identified with the intrinsic magnetic mo- 
ment of the wave packe1*i2iiLi&, is given by the expres- 
sion: 

M = ±|5R( M |vx (iJL_R)| tt ) = 
= ±^(r,*(u i \yx(i^-R)\u j )r lj ), 

where the sign is negative forelectrons and positive for 
holes. The operator v = y-^f- is the velocity operator 
corresponding to the degenerate part of the Hamiltonian. 

Written explicitly in component form and restricting 
our attention to holes, the magnetic moment is: 

M a = -I^NTMtf/s. (*^- ~ R 7 )}\^H- (42) 

e "/37 represents the antisymmetric tensor. It is straight- 
forward to prove that 



J7,*(u»|v|«j)»7j = I7 f vry 



1 de 
ftdk' 



(43) 



in which vy = (uj|v|uj). Therefore the second term in M 
is: 

M 2 = |^ V7?xR = _^ xR . (44 ) 
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The first term is 



in out 



Mi = -km xK ~l flt ^^ T,iv " !4r,! 



where 'out' means the sum runs over all bands outside 
the degenerate subspace, that is 1 7^ i,j. The first contri- 
bution exactly cancels M 2 , so the final result is 



where Hmtt has been defined in and g is the Lande 
g-factor. The first part of the Hamiltonian is Hd while 
the Zeeman term is H n . The Zeeman interaction between 
the spin and the magnetic field does not contribute to the 
velocity operator and therefore it does not contribute to 
the magnetic moment. The light-hole intrinsic magnetic 
moment in the spherical four-band model is given by the 
following expression: 



(46) 



Thus the magnetic moment can be expressed purely in 
terms of matrix elements connecting the degenerate sub- 
space to bands outside the subspace. 

We take as an example once again the light-hole mani- 
fold of the four-band Luttinger model in the spherical ap- 
proximation in the presence of a constant uniform mag- 
netic field. The Hamiltonian in this case is: 



H = H 



Lutt 



m 



(47) 



M = 



3efry 2 k 3 



(48) 



The magnetic moment is proportional to the expecta- 
tion value of the third component of the pseudo-spin and 
therefore to the helicity, as shown in (|39|l . Depending on 
the weight of each band in the wave packet the intrinsic 
magnetic moment can be positive or negative and if the 
bands are equally represented it will be zero. 
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